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1. Introduction 



We) = | V\ ) COS 6>i2 + 1 1/2) sin #12, (1.1) 

where #12 is the mixing angle. A more complete description would require the consideration 
of three neutrino generation mixing, but the available data do not allow to determine all 
the corresponding mixing parameters. 

If the neutrino mass states have also a non-vanishing magnetic moments, radiative 
decays 

v% — > Vj +7 (1-2) 

with rrii > m,- could be possible, as initially hypothesized in ||; the first searches for such 
decays were based on astrophysical considerations (see eg. || ) . The status of the decaying 
theory and phenomenology was summarized in (TO]. 
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Abstract: We present the results of the measurements performed in the occasion of the 
2001 total solar eclipse, looking for visible photons emitted through a possible radiative 
decay of solar neutrinos. We establish lower limits for the V2 and ^3 proper lifetimes above 
■^j" . 10 3 s/eV, for neutrino masses larger than 10~ 2 eV. 
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O ■ The data accumulated in the last few years in favor of neutrino oscillations (both solar, 

UJ, ^, [| and atmospheric |[], |4], |[) demonstrate that neutrinos have non vanishing masses 
(for recent reviews see [||, 0]). The most common interpretation of the oscillation data 
is based on two-generation mixing scenarios; the solar v e neutrinos are supposed to be a 
mixing of two mass eigenstates, v\ and ^2, with 7712 > m± (in the "normal hierarchy"): 
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The astrophysical neutrino lifetime lower limits are usually large (e.g. To/m > 2.8 x 10 15 
s eV _1 where tq is the lower proper lifetime limit for a neutrino of mass m, [pf|), but they 
are indirect and rather speculative limits. 

Much lower "semi-indirect" limits were deduced from the re-interpretation of solar and 
atmospheric neutrino data. Earlier attempts to explain the solar neutrino or atmospheric 
neutrino anomalies only in terms of neutrino decay have been ruled out by the existing 
evidence [12]; the present accepted explanations are based on neutrino oscillations, but do 
not exclude the hypothesis of neutrino decays. As an example, from the SNO data || a 
proper lower limit of To/m > 8.7 x 10~ 5 s eV _1 was deduced |l3|]. By analyzing all available 
solar neutrino data, other limits were obtained: r /m > 2.27 x 1CT 5 s eV" 1 for the MSW 
solution, and To/m > 2.78 x 10~ 5 s eV _1 for the vacuum oscillation solution of the solar 
neutrino problem (SNP) |14|j , or, following a different approach, r /m > 10~ 4 s eV _1 flip]] . 

Direct searches for radiative neutrino decays have been also performed. As an example 
we quote here the search for decay photons in the visible spectrum performed in the vicinity 
of a nuclear reactor Jl^j , yielding tq/ui lower limits in the range 10 -8 to nearly 0.1 s eV _1 , 
assuming neutrino relative mass differences Am/m between 10 -7 and 0.1. Recently, a 
search for 7 photons, using the Prototype Borexino Detector at Gran Sasso [|l7| reported 
To/m lower limits of 1.5 x 10 3 s eV _1 (assuming a polarization parameter a = —1 for the 
parent neutrino), 4.4 x 10 3 s eV _1 (for a = 0) and 9.7 x 10 3 s eV _1 (for a = +1). 

Total solar eclipses (TSO) represent a particular opportunity to look for radiative 
decays of solar neutrinos in the visible spectrum, during their flight from the Moon to the 
Earth, inside the shadow cone produced by the Moon. The first experiment based on this 
idea was performed in October 24, 1995 ]l8| , and a lower limit for the V2 proper lifetime To of 
about 10 2 s was obtained, assuming neutrino masses of few eV and Am^ = m\—m\ ~ 10 -5 
eV 2 . 

Some of us intended to perform measurements along this line during the 1999 TSE, 
in Romania: two experiments were prepared, one airborne and one at mountain altitude, 
but the weather conditions made the observations impossible |jl9|| . We could only Analise 
a video film recorded by a local television, obtaining V2 lower lifetime limits 1.8 x 10 -2 < 



r /m < 14.5 s eV" 1 gg, |21 



In this paper we present the results obtained from the analysis of our 2001 TSE ob- 
servations. 



2. Experimental data 

The July 21, 2001 TSE was observed from a location near Lusaka (Zambia) (14°56' lat. 
S, 28°14' long. E, and 1200 m a.s.l.), at about 8 km from the line of centrality. We used 
two instruments: a digital videocamera with an optical zoom 10 x and an additional 2x 
lens (which we will refer in the folowing as device "A" ) and a small Matsukov - Cassegrain 
telescope {<p = 90 mm, / = 1250 mm) coupled to a digital camera (device "B"). The fields 
covered by a single pixel where about 10" x 10" for images A and about 1.14" x 1.14" for 
B. The experimental data consist in 4149 frames obtained with the videocamera (data set 
"A") and 10 digital pictures obtained with the telescope (data set "B"). Fig. || shows two 
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of those images. For our analysis we summed the three color channels (Red, Green and 
Blue) of the images, thus recomposing the "white" light signal. The totality phase of the 
TSE in our observation site was about 3.5 minutes, so the displacement of the Sun behind 
the Moon is not negligible. As the signal of solar neutrino decays should be correlated with 
the direction to the center of the Sun, we calculated, for each frame in both data sets, the 
relative position of the Sun and, after determining by fit the center of the Moon disk, we 
computed the shift to be considered when summing the images in order to have the center 
of the Sun in the same pixel. 

Both instruments were calibrated at the Catania and Bologna Astronomical Observa- 
tories. In the exposure conditions of the eclipse, the number of visible photons required 
for producing 1 ADU (Acquisition Digital Unit) was 7.3 10 4 for instrument A and 8.9 10 2 
for device B. 



3. The simulation 

In order to extract physical information concerning a possible neutrino radiative decays 
from TSE data, a previous knowledge of the characteristics of the expected signal in manda- 
tory. For the analysis of the 2001 data we developed a full 3-dimensional Monte Carlo (MC) 



simulation [23|, based on the predictions of a recent version of the Standard Solar Model 
|24| ]. We first randomly chose a solar neutrino production reaction, and, consequently, a 
neutrino energy and the point of its creation inside the core of the Sun. Since we are inter- 
ested only in neutrino decays that may produce signals in our detectors, we then generate a 
random photon arrival direction, inside the angular acceptance of our devices, and a decay 
point, uniformly distributed along the photon path, between the observation point (the 
Earth) and the Moon. 

This procedure does not respect the kinematic probabilities of the simulated decay, so 
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we attribute to each MC event a weight according the decay angular probability density: 
dT mi 



oc 



dcos9 * (Am 2 



ij 



(m 2 + m 2 + mirrij) (1 + a cos 9*) . (3.1) 



In Eq. 3A rrii and rrij are the masses of the parent and respectively daughter neutrino (see 
Eq. L2), Arri'j = m 2 — m 2 , and 9* is the angle between the photon momentum and the spin 
of the initial neutrino, in its center of mass reference frame. The polarization parameter a 
varies from -1 (left-handed) and 1 (right handed) for Dirac neutrinos, and is for Majorana 
neutrinos. The kinematic weight for each MC event is obtained by integrating Eq. 3.1 over 
all photon directions that would lead to a signal in the same pixel, thus considering the 
experimental angular resolution. Note that such weights are dependent on the simulated 
device. 

This MC incorporates both the realistic geometry of solar neutrino production and 
decay, (as in ]25]j , where considering monoenergetic solar neutrinos the authors obtained 
an analytical prediction) and the standard energy spectrum predicted by the SSM, (as in 
our previous code |]l9| , pof where the neutrino source was approximated as pointlike). We 
assumed that m\ < mi < m% where mi, mi, are the masses of the v\, vi and ^3 mass 
eigenstates, respectively. We restrict our analysis to a two generation mixing scenario, 
assuming the present mass differences obtained from solar neutrino experiments, the LMA 
solution with Am 2 2 = 6 x 10 -5 eV 2 . Since SNO suggests also the presence of V3 in the solar 
neutrino flux, we considered also the mass difference measured by atmospheric neutrino 
experiments: Am 2 3 ~ Am^ = 2.5 x 10 -3 eV 2 . 

Fig. H shows the expected "luminosity curves" (the average luminosity versus the 
angular distance from the center of the Sun) from the simulation of (a) vi — > v\ + 7 and 
(b) ^3 — > v\ + 7 decays. The weights in Eq. |3.1| where calculated for different m\ values. 
The histograms in Fig. ||a suggest that for all neutrino masses, the expected signal is 
concentrated at small 9e angles (about 50 arcsec). The widths and shapes of the signals 
are sensitive to the mass assumed: the larger the mass, the narrower the signal band. In 
the case of Fig. ^b, the signal is broader (about 250 arcsec) and is less sensitive to the 
mass choice. 

One of the main results of the MC simulation consists in the determination of the 



global probabilities P of a solar neutrino decay according to Eq. 1.2, during its flight from 
the Moon to the Earth, and, produces a visible photon that reaches the detector. Those 
probabilities are both neutrino mass and instrument dependent. Consequently, assuming 
that an experiment detects iV 7 photons produced by neutrino radiative decays, the lifetime 
of the neutrino can be calculated from 

N, = P^S M t ob s (l - e"^) e~ tJ ^ , (3.2) 

where P are the probabilities estimated by the MC simulation, <I>j = & u s'm 2 On, (<fr u is 
the flux of solar neutrinos at the Earth (or Moon) and On the mixing angle) is the local 
flux of solar V{ mass eigenstate neutrinos, Sm is the area of the Moon surface covered by 
the analysis and t b s is the time of observation. {Pme) is the average time spent by solar 
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Figure 2: The expected shapes of the visible signals produced by the hypothesized solar neutrino 
radiative decay, assuming mi — 0.001 eV (solid histograms), 0.01 eV (dashed) histograms) and 
0.1 eV (dotted histograms). The squared mass difference is assumed to be 6 x 10~ 5 eV 2 (a) and 
2.5 x 10~ 3 eV 2 (b). In all cases a = -1. 

neutrinos inside the observation cone (one third of the flight time from the Moon to the 
Earth), and tsM is the time of flight of the neutrinos from the Sun to the Moon. 
A complete discussion of the MC and of its results may be found in |^] . 

4. Data analysis 

From the simulation described in the previous Section it follows that the expected visible 
signal from solar neutrino decays would have specific angular scales. A proper tool for 
investigating different scales of the eclipse images is the wavelet analysis. This technique 
gradually removes the contributions from various background sources as the diffraction of 
the coronal light on the borders of the Moon, the diffuse sky light, the ashen light (light 
reflected by the Earth on the surface of the Moon), etc. We used the simple Haar wavelet 
basis pq| . The n-order term of the decomposition is obtained by dividing the N xN pixels 2 
image in square fields of N/2 n x N/2 n pixels 2 and averaging the luminosity in each field; 
the averages are then removed and the resulting image, the n-order residual, can be used to 
obtain the (n+ l)-order term. Thus, each decomposition term results in an image in which 
objects of the corresponding scale are dominant, while the residuals contain information 
for smaller dimension scales. 

As the wavelet analysis requires a dyadic dimension of the field (the number of pixels 
on each border of the image is a power of 2), we retained from our images 64 x 64 pixels 
(from data A) and 512 x 512 pixels (from data B) defined around the position of the center 
of the Sun, and summed them in order to increase the signal to background ratios. 
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Figure 3: The luminosity profiles obtained from the raw data: a) data set A, b) data set B. 



The decay signal is searched for by averaging the luminosity of the images over "rings" 
centered on the position of the center of the Sun. There is no "central pixel"; so we 
have considered each of the four pixels adjacent to the image center as "central" and 
then averaged the obtained luminosity profiles. Such profiles, obtained from the raw data 
(before the wavelet decomposition) are shown in Fig. |3[ The difference between the two 
total images is due to the different CCD sensitivities, spatial resolution and optical features 
of the instruments. Data set A presents no clear structure, data B might contain some at 
relatively large 8. In order to check that the shape of the luminosity curve in data set B 
could be produced by the ashen light, we aligned the 10 digital images along the direction 
of the center of the Moon, and made a similar analysis on a full Moon pictured obtained 
with the same instrument. Figs. |] show this comparison. In order to enhance the contrast 
of the images, from both pictures we removed the average luminosity, thus obtaining the 
O^-order residual of the wavelet decomposition. 

The structures in Fig. ||a are similar to those in Fig^b; we should take into considera- 
tion that the Earth reflects the light of the Sun as a convex mirror, thus the central part of 
the Moon receives more light from the Earth than the rest of it (the relative excess in the 
raw TSE data is only about 2%). Instead the Sun illuminates the Moon uniformly. This 
observation suggests that we cannot simply remove the image of the full Moon from the 
data, as it would create a fake signal in the central part of the resulting image. Thus, one 
should develop a reliable model of the ashen light; alternatively one should use the wavelet 
decomposition. Note that the exposure conditions for the image of the full Moon were 
different than those during the eclipse, so the ADU values are not directly comparable. As 
we cannot determine which is the real contribution of the ashen light in set B, we can use 
the results only to determine a lower limit for the 1/3 lifetime. 
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Figure 4: Light luminosity profiles after removing the average luminosity of (a) the sum of images 
in data set B aligned with respect to the center of the Moon, and (b): an image of the full Moon 
obtained with the same instrument. 



(b) 




Figure 5: "White" light luminosity distributions of the fourth order wavelet term of (a) the 
summed images A and (b) B (centered on the Sun). 

4.1 Search for the V2 — ► v\ + 7 signal 

The expected signal from a U2 — > v\ + 7 decay, considering its MC estimated width p3| , 
should be better seen in the fourth order wavelet term of data sets A and B, as it corresponds 
to structures with about 40" - 60" width. Figs. |H| show the luminosity distributions for 
this term. Note that each bin is an average over 4x4 pixels in the case of data set A (Fig. 
5a), and over 32 x 32 pixels for set B (Fig. 5b). No central maximum is present in both 
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Figure 6: The 95% CL lower limits for the V2 proper lifetime, as function of the m Vl , obtained 
from data sets (a) A and (b) B. The results are valid in a neutrino mixing scenario with only 
two generations, and Ara^ i = 6 X 10 -5 eV 2 . The discontinuities in the proper lifetime limits for 
righthanded neutrinos originate in the MC probabilities and reflect the changements in the initial 
neutrino energy imposed by the condition of obtaining visible decay photons pointing to the Earth. 
Other relevant limits are also indicated (see text). 



data sets, and we can use these wavelet terms to determine lower lifetime limits for the 
investigated decay. 

The 95% CL lower limits for the v<± lifetime, in our case of no signal, are obtained by the 
substituting in Eq. 3^2 iV 7 with 3a of the forth order wavelet terms decomposition of the 
data, and considering sin 2 9ij = sin 2 612 — 0.74 (the LMA solution of the "Solar Neutrino 
Problem", jl], ||, ||]). They are shown with thicker lines in Figs. |6|a (data A) and |6]b (data 
B), assuming that V2 is a Dirac (lefthanded or righthanded) or a Majorana neutrino. The 
recent limits obtained from the Borexino Counting Test Facility [17] are also shown, for 
comparison. The arrows labelled "SNO" and "WMAP" indicate the lower neutrino mass 
limit reported by SNO [Q, [|, and the upper mass limit obtained by WMAP [27]. The limit 
obtained by the first TSE experiment [18| is indicated by the horizontal arrow; note that 
this limit, obtained using different physical hypotheses, is valid for neutrino masses of few 
eV. 

Neutrino lifetime values larger than our lower limits are not in conflict with the oscil- 
lation explanation of the solar neutrino deficit. The neutrino time of flight from the Sun 
to the Earth is about 500 s (in the laboratory frame of reference). The Lorentz boost for 
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Figure 7: The 95% CL lower limits for the proper lifetime, as function of the v\ mass, ob- 
tained from (a) set A (a) and (b) B. The solid, dot-dashed and dashed lines correspond to three 
different neutrino polarizations, a = 1, and -1, respectively. These results are obtained assuming 
Am 3,(2,i) = 2 - 5 x 10_3 ey2 and sin 2 6»i 3 ~ 0.1 

a solar neutrino with a mass of about 0.02 eV is 7 ~ 1.5 x 10 7 , so the fraction of ^2 that 
would decay into v\ + 7, assuming to ~ 60s (in the cm.) would be only ~ 5 x 10~ 7 . 

4.2 Limits on the v% — > v\p + 7 lifetimes 

As already shown, the search for v% — > v\ i + 7 signals is more difficult, as the ashen light 
could create a fake signature. Furthermore, the expected angular width of the expected 
signal is larger, so the wavelet decomposition could erase it. We still can compute 95% CL 



lower limits for the corresponding lifetimes, substituting A^ 7 in Eq. 3.2 with 3o"at 7 of the 
raw data A and B. The mixing angle #13 is not known, but it should be small; we assume 
sin 2 #i3 ~ 0.1. The lifetime limits obtained in those conditions are shown in Fig. [?]. 

5. Conclusions 

We analyzed two sets of digital images obtained during the June 21 st 2001 total solar 
eclipse, in Zambia, looking for possible radiative decays of solar neutrinos, yielding visible 
photons. 

Data set A consists in a large number of frames recorded with a digital videocamera; 
it has a relatively large integration time, but a modest space resolution. 
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Set B consists of 10 pictures taken with a digital camera coupled to a small telescope. 
Its time coverage is poorer than for set A, but it has a better space resolution and the 
instrument sensitivity was an order of magnitude better. 

The proper lower lifetime limits (95% CL) obtained for the vi — > v\ + 7 decays of 
lefthanded neutrinos range from To/m,2 — 10 s eV~ x to ~ 10 9 s eV -1 , for 10 -3 eV < m Ul < 
0.1 eV, see Fig. pi These limits are among the best obtained from direct measurements, 
demonstrating the potentiality of neutrino decay experiments during total solar eclipses 
(or possibly made in space, using the Earth as light absorber |2Ej). The lab. lifetime limits 
are about 10 7 times larger, thus the fraction of neutrino decays from the Sun to the Earth 
would be negligible. 

A similar analysis was made for a possible 1/3 — > ^2,1 +7 decay, assuming sin 2 6*31 ~ 0.1 
(the value of this mixing angle is not known). No signal compatible with a possible 1/3 — ► 
u 2,i + 7 is seen. The obtained 95% C.L. z/3 proper lifetime lower limits, for mi > 10 -2 eV 
and for a = —1, 0, are about two orders of magnitude lower than for the ^2, Fig- 

New observations, in better technical conditions, during forthcoming TSE's should be 
considered. 

An attempt along these lines was made during the December 2002 eclipse, but the 
weather conditions in South Africa did not allow any observation. We intended to use 
three portable telescopes, equipped with astronomy type CCD's. The sensitivity would 
have been about two orders of magnitude better than what reported in this paper. 
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